Symmetric achromatic low-beta collider interaction region design concept 
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We present a new symmetry-based concept for an achromatic low-beta collider interaction region 
design. A specially-designed symmetric Chromaticity Compensation Block (CCB) induces an angle 
spread in the passing beam such that it cancels the chromatic kick of the final focusing quadrupoles. 
Two such CCBs placed symmetrically around an interaction point allow simultaneous compensation 
of the lst-order chromaticities and chromatic beam smear at the IP without inducing significant 
2nd-order aberrations to the particle trajectory. We first develop an analytic description of this 
approach and explicitly formulate 2nd-order aberration compensation conditions at the interaction 
point. The concept is next applied to develop an interaction region design for the ion collider ring of 
an electron-ion collider. We numerically evaluate performance of the design in terms of momentum 
acceptance and dynamic aperture. The advantages of the new concept are illustrated by comparing 
it to the conventional distributed-sextupole chromaticity compensation scheme. 

PACS numbers: 29.20.db, 29.27.Bd, 41.75.-i, 41.85.Gy 
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I. INTRODUCTION 

In order to achieve the highest possible luminosity in 
a collider [l|, Q , the colliding beams must be focused to 
a small spot at the Interaction Point (IP). This tight fo- 
cusing is unavoidably accompanied by a large transverse 
beam extension before the beam enters the Final Focus- 
ing Block (FFB). The size of the required beam extension 
is determined by the desired degree of beam squeezing at 
the IP and the focal length of the FFB, which is closely 
related to the space between the IP and the nearest focus- 
ing quadrupole. The larger this distance, the greater the 
beam extension must be. Since the FFB's focal length 
depends on the particle's momentum, the large beam size 
inside the FFB leads to a large correlation between the 
particle's phase advance and its momentum Q. 

The problem with such a correlation is two-fold. First 
of all, it induces a large chromatic betatron tune spread. 
Since the available betatron tune space is limited by the 
beam resonances, the chromatic betatron tune depen- 
dence limits the ring's momentum aperture. Secondly, 
the chromatic dependence of the focal length causes chro- 
matic beam smear at the IP, which can even dominate 
over the beam size due to the emittance, significantly in- 
creasing the beam spot size at the IP and resulting in 
luminosity loss. Thus, an Interaction Region (IR) de- 
sign must employ sextupole magnets to compensate the 
chromatic effects However, the nonlinear sextupole 

fields generate 2nd- and higher-order aberrations to the 
particle trajectory at the IP, introduce nonlinear phase 
advance and limit the ring's dynamic aperture. Com- 
pensation of these nonlinear effects is one of the main 
challenges of an IR and, even more generally, collider 
ring design. A commonly used chromaticity compensa- 
tion technique is to install same-strength sextupoles in 
pairs with —I (minus the identity) transformation be- 
tween them to cancel their nonlinear kicks 0, HI ■ How- 
ever, this approach does not treat all of the sextupole- 



induced 2nd-order effects. Besides, organizing —I trans- 
formations for the sufficient number of sextupole pairs for 
simultaneous compensation of the horizontal and vertical 
chromaticities is difficult. The concept described below 
not only includes all of the —I transformation benefits 
but goes further by compensating all of the 2nd-order ef- 
fects together with the 2nd-order beam smear at the IP. 
It also offers a more compact compensation scheme. 

The IR design approach presented in this paper in- 
volves installation of a dedicated Chromaticity Compen- 
sation Block (CCB) next to the FFB In the 
CCB, certain symmetries of the beam orbital motion 
and dispersion are created using a symmetric arrange- 
ment of dipoles and quadrupoles. Two such CCBs are 
placed symmetrically around the IP. The symmetries of 
the beam orbital motion and dispersion combined with a 
symmetric quadrupole and sextupole arrangement in the 
CCBs then allow simultaneous compensation of the lst- 
order chromaticities and chromatic beam smear at the IP 
without inducing significant 2nd-order aberrations and 
therefore helping preserve the ring's dynamic aperture. 

A schematic IR layout is shown in Fig. [TJ The Beam 
Extension Section (BES) first expands the beam from 
its regular size in the arcs to the size required for fi- 
nal focusing. The beam next passes through the CCB, 
which creates in it an angle spread negatively correlated 
with the chromatic kick of the FFB, so that the FFB's 
chromatic effect is cancelled. The FFB then focuses the 
beam to the required spot size at the IP. The IR design 
is mirror-symmetric with respect to the IP. 



II. ANALYTIC DESCRIPTION 

In general, compensation of chromatic effects of a fo- 
cusing lattice requires introduction of dipole and sex- 
tupole magnets along the lattice. The dipole field 
changes the beam reference orbit and creates disper- 
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FIG. 1: Schematic IR layout consisting of Beam Extension 
Section (BES), Chromaticity Compensation Block (CCB), 
and Final Focusing Block (FFB). 



sion, i.e. orbit dependence on the particle energy. The 
sextupole field then introduces an additional energy- 
dependent focusing strength to counter the original chro- 
matic effect of the quadrupoles. 

The particle motion in the horizontal x and vertical 
y planes about the reference orbit is described by the 
following equations [Hj]: 



x" + {K 2 -n)x = K8 + f x , 
y" + ny = f y , 



(1) 



where the prime ' denotes a derivative with respect to 
the longitudinal coordinate s, K = 1/p is the curvature 
of the bent reference orbit, n = (dB y /dx)/(p/\e\) is the 
normalized quadrupole strength, and 6 = Ap/p is the 
particle's relative momentum deviation from the refer- 
ence value. The force terms f x and f y in Eq. ([T]) up to 
the 2nd order in S, x, and y are given by 



f x = -nSx + n s (x 2 - y 2 ), 
J y = n5y - 2n s xy, 



(2) 



where n s = (d 2 B y /dx 2 )/(p/\e\) is the normalized sex- 
tupole strength. The quantities K, n, n s , and, therefore, 
f x and f y are functions of s. 

The particle trajectory in the horizontal plane can be 
represented as a superposition of the dispersive compo- 
nent D 5 and the betatron component x b : 



x = DS + x b . 



(3) 



where the dispersion function D = dx/88 satisfies the 
equation 



D" + (K 2 - n)D = K. 



(4) 



Then the evolution of the betatron component is deter- 
mined by the following equations: 



x'l + (K 2 - n)x b = f x = -(n - 2Dn s )Sx b 

+ D(Dn s ~ n)S 2 1 - ' - 2 
y" + ny = f y = (n - 



±n s (xl -y~ 
2Dn s )Sy - 2n s x b y. 



(5) 



Sextupoles are usually placed in the ring arcs to com- 
pensate the energy dependence of the global betatron 
tunes. In that case, the chromaticity compensating sex- 
tupoles need to be configured to satisfy only two condi- 
tions associated with the Sx b and 5y terms on the right 
hand side of Eq. ([5]) , i.e. the integral effect of these terms 
over one revolution must be zero. The influence of the 
other three terms in Eq. ([5]) on the particle amplitude 
and global phase advance is zero when averaged over 
many revolutions. However, local chromatic effects in 
short sections or even at specific points like the IP where 
the beam experiences a strong transverse compression are 
generally not compensated on average in one turn. Com- 
pensation of these effects is a much more complicated 
problem, which requires taking into account all terms in 
Eq. ([S]). As one can see from these equations, the sex- 
tupole field component and dispersion give rise to the 
undesirable quadratic aberrations: 



Ax" oc n s [(DS) 2 + x\ 



- DnS 2 



Ay" 



-2n s x b y. 



(6) 



Compensating the impact of these terms on the beam 
size at the IP, where the beam size is supposed to be ex- 
tremely small, requires satisfying five different conditions 
in total. In order to explicitly state these conditions, we 
use the analytic perturbation approach below. 
We represent the betatron motion as 



x b = Xl + X, 

y = vi + y- 



(7) 



where the unperturbed linear parts x\ and y\ satisfy the 
homogeneous equations: 

x'( + (K 2 - n)xi = 0, 



Vi + nyi = 0, 



(8) 



while the perturbed trajectory components are given by 
x" + (K 2 - n)x - fx, 
y +ny = f y . 

We next express X\ and y\ as linear combinations of 
the two linearly- independent solutions u(s) and v(s) of 
Eq. ©: 

xi = a x u x (s) + b x v x (s), 

J/l = ClyUy(s) + byVy(s). 

We then use the parameter variation technique to solve 
for the perturbed part 

£ = a x (s)u x (s) + b x (s)vx(s), 
y = a y (s)u y (s) + b y (s)v y (s). 



(10) 



(11) 



Using the Cramer's rule, the formal solution for the per- 
turbed part can be written as 



X(s) = -^J-[v x [ Uxfxds-Ux [ Vxfxds], 

yvx Jo Jo 

y^ = 4r^ v v[ u vfv ds - u v I v v fyds], 

vv y Jo Jo 



(12) 
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where W x = u x v' x — u' x v x and W y are the Wronskians of 
the x and y basis functions u and v. 

Suppose that u and v are chosen such that at the IP 
u* = and v'* = 0, then the perturbed trajectory com- 
ponents at the IP are given by 



IP 



ttK [ u x f x ds = — 3_ f Ux f x ds, 
y x Jo u x J Q 



w. 

r = W y V vj UvfvdS 



(13) 



— , UyfydS. 

u y JO 



By performing the integration in Eq. (|13p over the whole 
beam path from the start of the beam extension after 
the arc to the IP, we obtain the perturbation of the par- 
ticle transverse position at the IP caused by the aber- 
rations from all of the involved bends, quadrupoles, and 
sextupoles. 

To find the next order correction to the unperturbed 
solution x\ and yi, in accordance with the iteration 
method, we plug the lst-order unperturbed solution x — 
D5 + x\ and y = y\ into the right-hand side of Eq. (|9|), 
obtaining the following equations for the 2nd-order cor- 
rection: 



nx2 = —(n — 2Dn s )5xi 
D(Dn s - n)5 2 + n s (x\ — yf), 
ny 2 = (n - 2Dn s )Sy 1 - 2n s x 1 y 1 , 



(14) 



where we neglected the term K 2 x on the left-hand side 
of the x equation, since the effect of the orbit curvature 
on beam focusing at high energy is very small compared 
to that of the quadrupolc field. Using Eq. (JT3J) , the 2nd- 
order correction to the particle's unperturbed trajectory 
at the IP is given by 



7j / u x [-(n - 2Dn s )5x 1 

u x JO 

+ D(Dn s -n)5 2 +n s {x 2 -yl)]ds, (15) 

7j / u y [(n - 2Dn s )5y 1 - 2n s x 1 y 1 ]ds. 

u v Jo 



2/2 = 



Consider an almost parallel beam after the beam ex- 
tension. The function u then describes the dominant 
("cos" -like) parallel component of the trajectory while 
the orthogonal solution v is associated with the small 
remaining angular spread ("sin" -like trajectory). Then, 
neglecting the angular divergence, one can use Eq. (|T0|) 
to approximate x\ and y\ in Eq. (|15j) as 



xi w a x u x {s), 
yi w a y u y (s), 



(16) 



x 2 



2/2 



—JZI u x [-(n - 2Dn s )Sa x u x 
< Jo 

D(Dn s ~ n)S 2 + n s (a 2 u 2 — a 2 u 2 )}ds, 
1 f* 

— / u y [(n - 2Dn s )Sa y u y 

u y JO 

^ifb g ^-^x ] * 



(17) 



Clearly, aberration compensation requires 



i* =0, 
2/ 2 * = 0. 



(18) 



Equation (jTTj) then leads to the following five compensa- 
tion conditions: 

/•# /• * 

(1) 2 / Dn s u 2 ds = / nu 2 ds, 
Jo Jo 



(2) 2 / Dn s u 2 ds 



:ds, 



(3) / D{Dn s - n)u x ds = 0, 



(19) 



(4) / n s u x ds = 0, 



Conditions (1) and (2) in Eq. (fT9|) are the usual condi- 
tions for compensation of the chromatic spread. Condi- 
tions (3)-(5) are the requirements for compensation of the 
2nd-order dispersion effect and sextupole beam smear. 

In the above analysis of the perturbed trajectory com- 
ponent, we so far ignored contribution of the angular 
spread, which is associated with the solution v. Even 
though this contribution is small compared to that of 
u, it may become dominant after compensation of the 
main aberration terms. Perhaps more importantly, the 
nonlinear phase advance associated with the angular 
spread may adversely affect the dynamic aperture. Us- 
ing Eqs. (fT0| and (fT5|) and neglecting the terms of the 
2nd order in v, contribution of the angular spread to the 
2nd-order perturbation at the IP can be written as 

~ \z I u x [-(n - 2Dn s )6b x v x 

u x Jo 

~\~ 2n s (d x u x b x v x (Xytty b y Vy)]ds, 



^2/2 ~ 1Z \ u v\-( n ~ 2Dn s )8byV y 

u v Jo 



(20) 



S Jo 

v x )]ds. 



obtaining for the 2nd-order trajectory correction at the 



Compensation of the aberration terms in Eq. (|20[) im- 
poses five more constraints in addition to those in 
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Eq. (HD: 



(1) 
(2) 
(3) 
(4) 
(5) 



Dn s u x v x ds = I nUxVxds, 
ii Jo 



2 / Dn s u y Vyds = 



n s u x v x ds = 0, 



n s u x u v Vyds = 0, 



n s u y v x ds = 0. 



HUyVydS, 



(21) 



Corrections beyond the 2nd order can be obtained by 
considering further iterations: 



(22) 



X — x 2 + x 3 + . . . , 

y = m + §3 + ■ ■ ■ 

Thus, the 3rd-order correction terms can be found by 
integrating the following equations: 

x'l — nx 3 n5 2 xi + (2Dn s — n)8x 2 

+ 2n s (xix 2 - 2/12/2) - \n"(D8 + xi) 3 
6 



+ n [(D5+x 1 ) 3 - 3(DS + x^yf], 

n5 2 yi - (2Dn s - n)Sy 2 
1 



y'l + ny 3 



(23) 



6 yi 



- 2n s (xiy 2 + 2/1X2) 
+ n [y 3 -3{DS + x 1 ) 2 y 1 ], 

where we introduced the octupole field component n = 
(d 3 B y /dx 3 )/(p/\e\). The higher-order terms should be 
studied to determine their effects on the beam smear at 
the IP and the ring's dynamic aperture. Equation (|23[) 
can be integrated similarly to Eq. (1141) : however, a com- 
plete analytic analysis clearly becomes rather cumber- 
some even in this 3rd-order case and one must employ 
a numerical approach. Nevertheless, one can still use 
Eq. (|23p to identify the symmetries helpful in compensat- 
ing the higher-order aberrations to guide the numerical 
studies. 



the five conditions of Eq. (fl9l) are reduced to just the 
first two, i.e. the linear chromaticity compensation con- 
ditions. While creating chromatic kick, the CCB does 
not generate the nonlinear aberrations associated with 
the 2nd-order effects of the dispersion and transverse 
beam sizes in the sextupole fields. All these terms corre- 
sponding to conditions (3)-(5) of Eq. (fH))) are automati- 
cally compensated inside the CCB due to the special op- 
tics design and lattice symmetry. Equation (p~9|) ignores 
the effect of the beam angular spread, which is admit- 
tedly small since the beam is assumed to be greatly ex- 
tended and almost parallel at the entrance into the CCB. 
However, this effect may adversely affect the ring's dy- 
namic aperture. The 2nd-order terms associated with the 
beam's angular spread are given by Eq. (|21[) . The sym- 
metry conditions in Eq. (|24p do not give any advantage 
in this case. The significance of these and higher-order 
terms needs to be examined in each specific case and 
any required compensation may be done by introducing 
additional sextupole and higher-order multipole families, 
perhaps, also with appropriate types of symmetry. 

We next test our symmetry concept by developing a 
conceptual IR design [TJ [l2| for the ion collider ring of 
the Medium-energy Electron-Ion Collider (MEIC) that 
is being proposed by Jefferson Lab [l[. Note that this is 
a proof-of-principle demonstration rather than a final IR 
design. It is used primarily to validate our symmetry- 
based chromaticity compensation concept. However, due 
to the modularity of our approach and the fact that the 
concept is independent of a particular FFB design, the IR 
design presented in this paper can be easily optimized to 
meet any specific nuclear physics detector requirements 
as discussed in 0, [lj| . 

Figure [2] shows a proposed layout of the MEIC com- 
plex [l|. The main MEIC parameters are summarized 
in Table HI As shown in Fig. [21 the electron and ion col- 
lider rings have geometrically-matching figure-8 shapes 
consisting of two 240° arcs connected by two straight 
sections crossing each other in the middle at 60°. There 
is a total of two IPs: one IP per straight. The rings are 
stacked vertically in the arcs but, at the IPs, the beams 



III. IR LINEAR OPTICS DESIGN 

Equation ()19|) is an analytic representation of the in- 
teraction region design requirements up to the 2nd order 
in effects on the particle trajectory. By imposing certain 
symmetries of the orbital motion, dispersion and mag- 
netic field components with respect to the center of the 
CCB, namely, 

u x (s) = pit x (-s), 

D(s) = ±D(-s), (24) 
n(s) = n(—s), 
n s (s) = ±n s (-s), 



Prebooster 




FIG. 2: Schematic layout of the proposed MEIC complex. 
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TABLE I: Main MEIC parameters. 



Parameter 


Unit 


Electrons 


Protons 


Momentum range 


GeV/c 


3-11 




20-100 


Optimization point 


GeV/c 


5 




60 


Polarization 


% 


~ 80 




> 70 


Collision frequency 


MHz 




748.5 




Particles per bunch 


10 10 


2.5 




0.42 


Beam current 


A 


3.0 




0.5 


Ap/p (rms) 


10~ 3 


0.7 




0.3 


Bunch length (rms) 


mm 


7.5 




10 


Norm, x emit, (rms) 




54 




0.35 


Norm, y emit, (rms) 


pLm 


11 




0.07 


A. at IP (K) 


cm 




10 




P y at IP (ft) 


cm 




2 




Collider ring length 


m 




1340 




Luminosity per IP 


cm _2 s _1 


5. 


6 x 10 


33 



e s a s q s q 

♦ A ♦ A ♦ 



FFB 




FIG. 3: A conceptual drawing of an IR with a CCB that 
has an odd-symmetry horizontal betatron trajectory Xb, an 
even-symmetry vertical betatron trajectory yb, and an even- 
symmetry horizontal dispersion D x . It shows a possible ar- 
rangement of dipoles, quadrupoles Q, and sextupoles S that 
satisfies the magnetic field symmetry requirements and pro- 
duces the desired symmetries of Xb, yb, and D x . 



cross in the horizontal plane. This is done by keeping 
the electron ring flat and bringing the ion beam to the 
electron beam's plane using vertical chicanes located at 
the ends of the arcs. 

Since the IR design concepts are similar for the two 
collider rings, below we focus our discussion on the ion 
collider ring, which is perhaps more challenging due to its 
larger natural linear chromaticities [l|. Two IRs are in- 
corporated into the ion ring's two straight sections. The 
IR sections are assumed identical and the ring is arranged 
in a two-fold symmetric way. The ring's horizontal v x and 
vertical v y betatron tunes are set to 23.273 and 21.285, 
respectively, by adjusting the phase advance outside of 
the IRs. The ring lattice used in the below simulations is 
fairly basic; it docs not include components such as RF 
cavities, detector solenoids, vertical chicanes, spin control 
devices, injection/extraction elements, etc. However, the 
details of the ring design outside of the IR sections have 
little impact on the results of this paper and their de- 
scription can be found in p|. Therefore, below we focus 
our attention on the IR design. 

A conceptual drawing of an IR with a CCB that has 
an even symmetry of the dispersion and an odd sym- 
metry of the horizontal betatron trajectory is shown in 
Fig. [3] If one then uses an even symmetry of the sex- 
tupole fields, conditions (3)-(5) of Eq. (fll?)) are automati- 
cally satisfied. The two remaining chromaticity compen- 
sation conditions are satisfied by adjusting the fields of, 
at least, two sextupole families. This can be attained us- 
ing the difference in the behavior of the horizontal and 
vertical /3-functions. 

In Fig. [3J there are two identical bends at the begin- 
ning and at the end of the CCB, which generate and 
then suppress the dispersion while symmetric quadrupole 
optics ensures the appropriate symmetries of the beta- 
tron and dispersive orbital components with respect to 
the center of the CCB. When designing the ion interac- 



tion region of an electron-ion collider, one has to keep in 
mind that it has to match the footprint of the electron 
interaction region. In the electron ring, since the CCB 
dipoles are located in regions with large /3-function val- 
ues, their maximum bends must be limited to avoid their 
cmittancc-dcgrading impact. In the ion ring, on the other 
hand, it is advantageous to have strong bends to produce 
a large dispersion thus reducing required sextupole fields 
and their nonlinear effects. These bend requirements are 
somewhat contradicting. Our solution is to use alternat- 
ing but still symmetric bends in the ion CCB. One can 
then generate a large dispersion while having the flexi- 
bility of the footprint geometry. 

Since the chromaticity correction scheme is indepen- 
dent of particular BES and FFB implementations, we 
first focus on the CCB design. Figure @] shows CCB 
optics developed for the MEIC ion collider ring fol- 
lowing the concept illustrated in Fig. [3] Note that, 
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FIG. 4: CCB optics. 
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FIG. 5: Complete IR optics of the MEIC ion collider ring 
The IP is located at the symmetry point in the middle. 



collider mode of operation. Since the injected beam's 
emittance is typically much greater than after acceler- 
ation and cooling, keeping the collider-mode optics at 
injection may result in an unacceptably large beam size. 
Therefore, to keep the magnet apertures reasonable, the 
maximum /3-functions occurring in the IR section may 
need to be reduced at injection and during acceleration 
and then restored for the collider operation. The de- 
scribed IR design allows a straightforward implementa- 
tion of such a "/3 squeeze" . Since entering and exiting the 
CCB the dispersion is suppressed and the a-functions are 
almost zero, the CCB optics is compatible with any ini- 
tial /3-functions without any adjustment required as long 
as a XtV have appropriate close-to-zero values. Thus, the 
size of the /3-functions inside the IR can be controlled by 
adjusting the /3-function values at the end of the BES. 
This requires only changing the optics of the BES with 
the rest of the ring optics intact. 



in comparison to the simplified schematic of Fig. |3l a 
larger number of quadrupoles were used to better con- 
trol the beta function behavior. The CCB is composed 
of eight symmetrically-arranged alternating bends with 
seven quadrupoles placed symmetrically between them. 
The sum of the absolute values of all dipole bending 
angles is 480 mrad, while, due to the alternation of 
their bending directions, the net CCB bending angle is 
120 mrad. The quadrupole strengths are adjusted to pro- 
duce a total CCB transfer matrix that meets the symme- 
try requirements of Eq. 

For the conceptual study presented in this paper, a 
simple quadrupole doublet is used for the FFB. As shown 
in Table HI the design values of the horizontal and vertical 
/3-functions at the IP (/3* ) are 10 and 2 cm, respectively. 
The distance from the IP to the front face of the near- 
est FFB quadrupole is assumed to be 7 m, which, in 
combinations with the /3* y values, determines the max- 
imum values of the /3-functions in the IR region. The 
FFB optics design also determines the Twiss function val- 
ues at the entrance into the CCB. The FFB quadrupole 
strengths, lengths and spacing are adjusted so that the /3- 
functions have roughly equal values at the entrance into 
the CCB and the a-functions {a x , v = — (d/3 x , y /ds)/2) 
are equal to zero in the middle of the CCB. A BES con- 
sisting of 7 quadrupoles then matches the CCB optics 
to the ring's regular lattice. This involves a substantial 
beam extension and therefore requires an appropriately 
large amount of longitudinal space. The IR is arranged 
symmetrically with respect to the IP. The complete IR 
optics of the MEIC ion collider ring is shown in Fig. [5] 
The detector solenoid and its coupling compensation el- 
ements are not included in this design. Note that the 
CCB bends on the opposite sides of the IP are reversed 
making the dispersion anti-symmetric with respect to the 
IP. This is also done for the purpose of electron and ion 
IR footprint matching. 

The optics shown in Figs. 0] and [S] corresponds to the 



IV. MOMENTUM ACCEPTANCE AND 
DYNAMIC APERTURE 

A. Chromaticity compensation and momentum 
acceptance 

In accordance with our chromaticity compensation 
concept, two sextupole pairs are inserted in each CCB. 
The sextupoles in each pair are identical and are placed 
symmetrically with respect to the center of the CCB. The 
sextupoles are shown in Figs. @] and [5] by the shorter bars. 
The sextupole positions are chosen at the points where 
the dispersion is near maximum and there is a large dif- 
ference between the horizontal and vertical /3-functions. 
The two parameters corresponding to the strengths of 
the sextupole pairs are adjusted to compensate the hor- 
izontal £ x and vertical t; y natural linear chromaticities. 

In this collider ring design [l[, the total values of £ x 
and £ y are -278 and -268, respectively. Contributions of 
the IRs to £ x and £ y are -254 (91.4%) and 245 (91.3%), 
respectively. Therefore, contribution of the remainder of 
the ring is almost negligible and can also be compensated 
locally by the CCBs. The two sextupole families are used 
to adjust the slopes of the chromatic horizontal and ver- 
tical betatron tune curves to zero. The chromatic tune 
dependence before and after the compensation is shown 
in Fig. [5] The horizontal and vertical tune variations 
are less than 0.02 and 0.03, respectively, within Ap/p of 
about ±4 x 10~ 3 . 

An effective method to evaluate the nonlinearity of a 
particle's dynamics is to obtain a frequency map [l5[ 
from particle tracking. At the design 60 GeV/c proton 
beam momentum of the MEIC ion collider ring, the max- 
imum horizontal and vertical rms beam sizes are 3.2 and 
1.6 mm, respectively, due to a relatively large detector 
space of ±7 m. This makes it more challenging to obtain 
a sufficiently large dynamic aperture in the horizontal di- 
mension. Therefore, the frequency map is first computed 
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FIG. 6: Momentum dependence of the horizontal v x (red) and 
vertical v y (blue) fractional betatron tunes before (solid lines) 
and after (dashed lines) the linear chromaticity compensation. 

in the [x — Ap/p) phase space as shown in Fig. [7] The 
map is obtained by tracking particles for 2000 turns us- 
ing a well-established accelerator simulation software El- 
egant [l6| . All ring components are modeled as canonical 
kick elements with exact Hamiltonians retaining all or- 
ders in momentum offset. The magnet fields are approx- 
imated as hard-edge and the lattice is assumed perfect, 
i.e. containing no alignment or field errors. The simu- 
lation does not take into account effects such as beam- 
beam, intra-beam scattering (IBS), etc. The particles are 
launched parallel to the beam axis at the entrance into 
one of the CCBs where the horizontal and vertical (3- 
functions are about 1330 and 1350 m, respectively. The 
color in Fig. [7] reflects the tune change in terms of the 
tune diffusion d defined as d = log(Ai^ + AVy) where 
Av XjV is the tune change from the first to the second half 
of the simulation in the horizontal and vertical planes, 
respectively. The diffusion index d is one of the best cri- 
teria [l3| to determine the particle's long term stability 



and enables one to see the nonlinear behavior, i.e. the 
more negative the diffusion, the smaller the nonlinear ef- 
fects are. 

Two conclusions can be drawn immediately from the 
frequency map in Fig. [7j 

1. the momentum acceptance can easily reach Ap/p of 
±4 x 10~ 3 with only the linear chromaticity com- 
pensation, 

2. the uniform color distribution means that there are 
no strong resonant perturbations in the particle 
tracking. 

The simulation was terminated at Ap/p = ±4 x 10 -3 , 
which is about HAas (an rms momentum spread as of 
about 3 x 10 -4 is expected after cooling) and demon- 
strates a sufficient momentum acceptance. Thus, with- 
out any compensation of higher-order chromaticities, the 
symmetry concept immediately results in an excellent 
momentum acceptance. Figure [7] also indicates that the 
horizontal dynamic aperture size is about ±10 mm, which 
is reasonable given the large compensated values of the 
natural chromaticities and the fact that there was no 
nonlinear optimization after the linear chromaticity com- 
pensation. However, due to the large beam extension re- 
quired to achieve the ambitiously small /3* values, clearly, 
further optimization of the dynamic aperture using mul- 
tiple sextupole and octupole families is required. One ap- 
proach to such optimization is presented in Section [IV CI 

Figure [8] shows the tune footprint corresponding to the 
same tracking as the frequency map in Fig. [7l As with 
the frequency map, the color of the tune footprint reflects 
the tune change in terms of the diffusion index. The lines 
indicate betatron resonances up to the 3rd order; higher- 
order resonances are not shown in this plot. Since the 
tunes in Fig. [8] are calculated for particles with initial 
coordinates in the (x — Ap/p) phase space, the vertical 
tune variation arises only from the chromatic tune depen- 
dence, which is around 0.03 within Ap/p of ±4 x 10~ 3 as 




FIG. 7: Frequency map in the (x — Ap/p) phase space. The 

color reflects the tune change^in terms of the tune diffusion d FIG. 8: Betatron tune footprint with the color indicating the 
defined as d = \og(Au x + Au y ). tune diffusion d. 
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shown in Fig. [6l The horizontal tune variation is about 
0.14, which is significantly greater than the chromatic 
tune change of 0.02 corresponding to Ap/p of ±4 x 10~ 3 
in Fig. [5] This indicates that, after the linear chromatic- 
ity compensation, the horizontal betatron motion pro- 
duces a substantial tune change, which can drive parti- 
cles close to or into a resonance mode and cause particle 
loss. The tune change due to the particle's transverse 
motion is discussed in detail in Section IIVCI Overall, 
the tune footprint immediately reveals the distribution 
of the particle tunes, which allows one to understand the 
tune trend due to the particle motion and subsequently 
adjust the design tunes to avoid approaching or crossing 
strong resonances. 

An important feature of our symmetry concept is com- 
pensation of the chromatic beam smear at the IP. It is 
attained together with the linear chromaticity correc- 
tion when two CCBs are placed symmetrically around 
the IP. One way to demonstrate this is by using the 
chromatic amplitude functions W XjV , which are directly 
related to the chromatic derivatives of the /3-functions 
df3 x ,y/dd fl8j |. Figure E] shows the W x>y functions plotted 
over the whole ion collider ring length before and after 
the linear chromaticity correction, respectively. The fact 
that W XyV are close to zero at the IP indicates that the 
chromatic beam smear is compensated. 



B. Comparison to the conventional 
distributed-sextupole scheme 

To demonstrate the advantages of the proposed sym- 
metry concept, we compared its performance with that of 
a more conventional "distributed-sextupole" (DS) chro- 
maticity compensation scheme based on introducing two 
sextupole families in the arc FODO cells. To make a fair 
comparison to the CCB scheme, the CCBs were taken 
out of the ring optics to eliminate their chromatic con- 
tribution, and the BESs were connected directly to the 
FFBs. Additional FODO cells and two more sextupole 
families were introduced in the regions between the arcs 
and BESs for higher-order chromaticity correction. The 
phase advances between the sextupoles and between the 
sextupoles and IP were adjusted to cancel the sextupole 
nonlinear kicks and to optimize the 2nd-order chromatic- 
ity compensation as discussed in Q . These modifications 
to the linear optics resulted, in particular, in higher be- 
tatron tunes but the initial horizontal £ x an d C-y vertical 
natural chromaticities of -199 and -262, respectively, were 
still smaller than in the CCB scheme. 

First, for the linear chromaticity compensation, two 
sextupole families were introduced in the arcs. Following 
the same guidelines as in the CCB case, the sextupoles 
were placed at locations with nearly maximum dispersion 
and a large difference between the horizontal and vertical 
/3-functions. With 52 FODO cells in the two arcs of the 
ion collider ring, one might expect that such a globally 
distributed sextupole scheme can effectively reduce the 
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FIG. 9: Horizontal and vertical chromatic amplitude func- 
tions W x ,y plotted over the whole ion collider ring length be- 
fore (top) and after (bottom) the linear chromaticity compen- 
sation. 



sextupole strengths with a consequent suppression of the 
nonlinear effects introduced by the sextupole fields. How- 
ever, the /3-functions of a regular arc FODO cell are much 
smaller than those in the CCBs. Hence, even though the 
DS scheme has a total of 104 sextupoles compared to only 
16 sextupoles in the CCB scheme, the sextupole strengths 
in the DS scheme are still much stronger than those in 
the CCB scheme. Table iHl summarizes the parameters of 
the two different linear chromaticity correction schemes 
including the values of the optical functions at the sex- 
tupole locations and the sextupole strengths. 

Figure [TU] compares the chromatic dependence of the 
fractional betatron tunes before and after the linear chro- 
maticity compensation using the two sextupole families 
in the arc FODO cells. The slopes of the horizontal 
v x and vertical v y chromatic betatron tune curves are 
zero at Ap/p = after the correction, which reduces the 
tune variation and somewhat extends the momentum ac- 
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TABLE II: Comparison of the optical functions at the sex- 
tupole locations and of the sextupole strengths d 2 B y /dx 2 for 
60 GeV/c protons in cases of the DS and CCB-based linear 
chromaticity correction schemes. 



Parameter 


DS scheme 


CCB scheme 






Sextupole family 






1 


2 


1 


2 


# of sextupoles 


52 


52 


8 


8 


A* (m) 


13.8 


6.0 


1614 


2.3 


I3 y (m) 


6.0 


13.8 


944 


1396 


D x (m) 


1.19 


0.8 


0.9 


1.3 


d 2 B y /dx 2 (T/m 2 ) 


4000 


-7000 


238 


-288 



ceptance. However, the second derivatives of the cor- 
rected chromatic curves (the 2nd-order chromaticities) 
start dominating the tune spread and still limit the mo- 
mentum acceptance. Therefore, we introduce two addi- 
tional families to a total of 128 sextupoles close to the 
IR to compensate the 2nd-order chromaticities. For the 
optimal compensation, we set the phase advance between 
the sextupoles and IP following the approach described 
in The strengths of the new sextupole families were 
comparable to those of the two original families in Ta- 
ble El 

Figure QT] shows the chromatic tune dependence after 
combining the linear and 2nd-order chromaticity com- 
pensations. Figure QT] indicates that the momentum ac- 
ceptance is extended further in comparison to Fig. 1101 
However, the strong sextupole fields give rise to even 
higher-order chromaticities still leading to an impracti- 
cally large tune spread. The frequency map and tune 
footprint for this 2nd-order DS chromaticity compensa- 
tion scheme are also calculated using particle tracking 
and plotted in Figs. [12] and [131 respectively. Because of 
such a large chromatic tune spread, the particles expe- 
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FIG. 11: Chromatic dependence of the fractional betatron 
tunes before and after the chromaticity compensation up to 
2nd order. 



rience numerous resonances, which limit the momentum 
acceptance. 



C. Dynamic aperture with the CCB scheme 

A ring's property that is equally important to its mo- 
mentum acceptance is its dynamic aperture (DA), i.e. 
a region in the transverse plane of stable particle mo- 
tion. Therefore, we study the ion collider ring's DA for 
the case of the CCB-based chromaticity compensation 
scheme. The DA is explored by tracking particles for 
1000 turns while progressively increasing their initial hor- 
izontal (x) and vertical (y) amplitudes until the boundary 
between survival and loss is found. Figure [14] shows the 
on-momentum particle's DA at the entrance into one of 
the CCBs. The DA reaches ~ 10 mm horizontally and 
~ 18 mm vertically, which is quite reasonable given the 
large natural chromaticities and the fact that there was 
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FIG. 10: Chromatic dependence of the fractional betatron 
tunes before and after the linear chromaticity compensation 
using the DS approach. 
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FIG. 12: Frequency map in the (x — Ap/p) phase space for 
the 2nd-order DS scheme with the color indicating the tune 
diffusion d. 
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FIG. 13: Tune footprint for the 2nd-order DS scheme 
the color indicating the tune diffusion d. 
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pendence is shown in Fig. 1151 v x is plotted as a function 
of x and y at the top while v y is plotted as a function of 
x and y at the bottom. 

The simulations show that, after the sextupole com- 
pensation, there is a large amplitude-dependent tune 
shift. This is due to the fact the nonlinearity of the sex- 
tupole fields causes stronger focusing of the particles with 
larger transverse amplitudes. As seen from Fig. 1151 v x 
and v y change from the design values by 0.08 (the red 
line in the top plot) and 0.04 (the blue line in the bot- 
tom plot), respectively, within x and y of about ±10 mm. 
Such a large tune deviation can easily drive particles to 
approach or cross betatron resonances, which inevitably 
causes their loss. The tune shift with amplitude limits 
the ring's DA, especially in the horizontal direction due 
to the large horizontal beam extension in the ion collider 
ring (the maximum horizontal rms beam size a x along the 
ring is 3.2 mm). Figure [151 also shows that there is very 
little coupling between the two transverse dimensions: 



no nonlinear optimization after the linear chromaticity 
correction. However, due to the large beam extension 
required to achieve the ambitiously small /3* values at 
the MEIC (see Table J)), these horizontal and vertical DA 
sizes correspond to only ~ 4ct x and ~ 15a y , respectively. 
Thus, optimization of the DA using multiple sextupole 
and octupole families is required. 

For the on-momentum particle, there is clearly no tune 
variation due to the chromatic tune dependence. The 
tune change comes from the nonlinearity of the trans- 
verse betatron motion caused by the nonlinear fields, in 
particular, of the sextupoles. Thus, we explore the de- 
pendence of the betatron tunes on the transverse motion. 
For a particle with initial transverse amplitudes x and y, 
its horizontal v x and vertical v y tunes are calculated by 
tracking it for a number of turns. The obtained tune de- 
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FIG. 14: Dynamic aperture in the (y—x) space at the entrance 
into one of the CCBs for the on-momentum particle. 
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FIG. 15: Dependence of the horizontal p x (top) and vertical v y 
(bottom) betatron tunes on the initial transverse amplitudes 
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v x is almost independent of the vertical motion (the blue 
line in the top plot of Fig. [T5]) while v y is unaffected by 
the horizontal motion (the red line in the bottom plot of 
Fig. US]). 

A straightforward mechanism to compensate the 
amplitude-dependent tune shift caused by the sextupoles 
is to introduce octupoles. Depending on the octupole lo- 
cations, both the amplitude dependent tune shift and 
higher-order chromaticities can be controlled by adjust- 
ing the octupole strengths. Since the momentum ac- 
ceptance is large enough after the linear chromaticity 
correction with the CCBs, we focus on compensation 
of the amplitude-dependent tune shift. Therefore, oc- 
tupole families should be placed in dispersion-free regions 
to leave the chromatic correction unaffected. Also, to 
reduce the required octupole strengths, they should be 
placed at large /3-function points. Following these consid- 
erations, three families of octupoles are introduced in the 
BES and FFB regions. It was the most convenient for us 
to find a solution for the octupole strengths that provide 
compensation of the lst-order amplitude dependent tune 
shifts, u sing the arbitrary-order expansion tools of COSY 
Infinity |19| . After this compensation, the horizontal and 
vertical tune variations are again obtained as functions 
of the initial transverse amplitudes x and y from particle 
tracking and are plotted in Fig. 1161 Compared to the tune 
changes of 0.08 and 0.04 within the transverse amplitude 
of about ±10 mm in Fig. [151 the horizontal and vertical 
tune changes are now 0.04 and 0.03, respectively, within 
an amplitude range of over ±15 mm. Consequently, the 
DA is increased as shown by the blue line in Fig. [T7] in 
comparison to the red line without the octupole opti- 
mization. 

After the octupole compensation, the limitation on the 
DA comes from the effect of the 2nd-order amplitude- 
dependent tune shift, which starts dominating the tune 
behavior in Fig. [TB] This provides a guideline that the 
DA can be further enhanced by considering and opti- 
mizing the 1st-, 2nd-, and even higher-order amplitude- 
dependent effects. One possibility for future studies is to 
use a genetic algorithm to optimize the DA directly by si- 
multaneously adjusting multiple sextupole and octupole 
families [2(| ■ Note that magnet errors are not included in 
the above simulations. Addressing the question of error 
tolerances is the subject of our future studies. 



V. CONCLUSIONS 

We proposed a new symmetry-based interaction re- 
gion concept for a high-luminosity collider ring design. 
Such a symmetric design allows simultaneous compen- 
sation of (a) the betatron tune spreads due to the lst- 
order chromaticities, (b) chromatic beam smear at the IP, 
i.e. chromatic spread of focal point, and (c) other 2nd- 
order aberrations to the particle trajectory at the IP and 
on average over one turn, such as the 2nd-order disper- 
sion effect and sextupole-generated beam smear due to 
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FIG. 16: v x (top) and v y (bottom) as functions of x and y 
after compensation of the lst-order amplitude dependent tune 
shifts. 



the beam size. According to this concept, two specially- 
designed symmetric Chromaticity Compensation Blocks 
(CCBs) are placed symmetrically around an IP. Each 
CCB induces an angle spread in the passing beam such 
that it cancels the chromatic kick of the final focusing 
quadrupoles. We developed an analytic description of 
this approach and explicitly formulated 2nd-order aber- 
ration compensation conditions at the IP. We verified our 
concept by developing a specific interaction region design 
and applying our chromaticity compensation scheme to 
a challenging case of a high-chromaticity small-/?* MEIC 
ion collider ring [lj. We then compared performance of 
our symmetric CCB approach with that of the conven- 
tional distributed-sextupole chromaticity compensation 
scheme. The CCB approach resulted in a much greater 
momentum acceptance with a smaller number of weaker 
sextupoles. The dynamic aperture after the chromaticity 
correction is reasonably large, especially given the large 
natural chromaticities. We identified the factors limiting 
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FIG. 17: Dynamic aperture in the (y— x) space at the entrance 
into one of the CCBs for the on-momentum particle without 
(red) and with (blue) the octupole compensation of the 1st- 
order amplitude-dependent tune shift. 



the dynamic aperture, namely, the 1st- and higher-order 
amplitude-dependent tune shifts, and demonstrated an 
approach to dynamic aperture optimization using oc- 
tupole families. Further optimization may be needed de- 
pending on the specific requirements. Evaluation of the 
error impact will be addressed in the future studies. 
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